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We have demonstrated that the middle-infrared (mid-IR) peaks in the optical conductivity 
spectra of CeX^ (X = Pd, Sn, In) can be explained by first-principle band structure calculation 
with the spin-orbit interaction. The mid-IR peak shapes in these materials are not identical 
to one another: CePd3, CeSn3, and Celn3 have a triple-peak structure, double-peak structure 
and broad single-peak structure, respectively. These peaks can be theoretically explained by 
the optical transition from the occupied state to the spin-orbit splitted Ce 4/ state. This result 
indicates that the mid-IR peaks originate from the simple band picture with the Ce 4/ state 
near the Fermi level, not from the conventional cf hybridization gap based on the periodic 
Anderson model. 
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Heavy fermion as well as mixed valent materials such 
as cerium (Ce), ytterbium (Yb), and uranium com- 
pounds commonly have a characteristic peak structure 
(namely, the mid-IR peak) in the middle-infrared (mid- 
IR) region below 1 eV in their optical conductivity [<r(o;)] 
spectra. 1 ) The origin of the mid-IR peak has remained a 
long debated issue so far. Recently, two proposals have 
been made in different studies regarding the origin of the 
mid-IR peak; one attributing it to the optical transition 
between the bonding and antibonding state of the hy- 
bridization gap between conduction and local 4/ elec- 
trons, namely cf hybridization gap, predicted by the 
periodic Anderson model (PAM), 2 ~ 5 ) and the other to 
a peculiar band structure in which the 4/ state is lo- 
cated near the Fermi level (£ F )- 6) The for mer study 
presented a persuasive scaling plot of the energy of the 
mid-IR peak versus the hybridization intensity. 5 -* Actu- 
ally, the cf hybridization plays an important role in the 
physics of heavy fermion systems. 7 ) The cf hybridiza- 
tion band predicted by PAM has been directly observed 
with Ce Ad — 4/ resonant angle-resolved photoemission 
spectroscopy of CeCoGei^Sio.s that is a typical heavy 
fermion compound. 8 ) On the other hand, the latter study 
pointed out that the mid-IR peak of YbALj can be ex- 
plained in terms of a simple band picture without ref- 
erence to the cf hybridization gap but they shifted the 
Yb 4/ level near the Fermi level. 6 ) 

The mid-IR peak is mainly observed in Ce and Yb 
compounds. However the shape of the mid-IR peak is 
different from one another. Yb compounds have a single- 
peak structure due to the optical transition from the Yb 
4/7/2 state because the transition from the Yb 4/ 5 / 2 
state does not appear in the mid-IR region but in the 
near-infrared region above 1 eV due to a spin-orbit split- 
ting (SOS) energy of about 1.4 eV. 9 ) On the other hand, 
almost all Ce compounds have a double-peak structure 
with an energy split of about 0.25 eV originating from 
the splitting of the Ce 4/5/2 an d 4/ 7 / 2 states due to 



the spin-orbit interaction (SOI). 10 - 1 For instance, CeSn3 
that is a typical mixed valent material has a double-peak 
structure. 11 -* In CePd3 that is another mixed valent ma- 
terial with the identical CusAu-type crystal structure, 
there is a triple-peak structure with equally spaced en- 
ergy splitting of 0.25 eV in the mid-IR region. 12,13 ) An- 
other isostructural material, Celn3, which is in the an- 
tiferromagnetic ground state with the Neel temperature 
(Tn) of about 10 K, also has a broad peak in the mid- 
IR region, not only in the paramagnetic phase above but 
also in the antiferromagnetic phase below TV 11 ) If the 
mid-IR peak of Celn3 originated from the cf hybridiza- 
tion gap, Celn3 should not have a mid-IR peak below Tn 
because the Ruderman-Kittel-Kasuya-Yoshida (RKKY) 
interaction predominates over the cf hybridization in the 
ground state. To clarify the origin of the different mid-IR 
peak shape in CeY3 (X — Pd, Sn, In) is important for 
the identification of the origin of the mid-IR peaks in Ce 
compounds. 

In this Letter, to investigate the different shapes of the 
mid-IR peaks of CeY3 as well as the origin of the mid-IR 
peaks, we calculate the a(ui) spectra of these materials 
from the first-principle band calculation with SOI. The 
result shows that the mid-IR peak can be explained in 
terms of the optical transition derived from the band 
calculation. Therefore, it is concluded that the mid-IR 
peaks originate from the characteristic band structure, 
not from the conventional cf hybridization gap. 

The band structure calculation was performed us- 
ing the full potential linearized augmented plane wave 
plus local orbital (LAPW+lo) method including SOI 
implemented in the Wien2k code. 14 ) CeY 3 (X = Pd, 
Sn, In) forms a CusAu-type cubic crystal structure 
(PmSm, No. 221) with lattice constants of 4.1280 A, 15 ) 
4.7423 A, 15 ) and 4.6876 A, 15 - 16 ) respectively. The non- 
overlapping muffin-tin (MT) sphere radii values of 2.50 
Bohr radius were used for both the Ce and X atoms in 
CeA^. The value of RuT^max (the smallest MT radius 
multiplied by the maximum k value in the expansion of 
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(a) With spin-orbit interaction 



(b) Without spin-orbit interaction 




Fig. 1. Calculated electronic structure of CePd3 with (a) and 
without (b) the spin-orbit interaction. The band with solid cir- 
cles in (a) is the initial state of the mid-IR peaks. 



plane waves in the basis set), which determines the ac- 
curacy of the basis set used, was set to 7.0. The total 
number of Brillouin zones was sampled with 40,000 k- 
points. The band structure of LaPd3 (lattice constant = 
4.235 A 17 ') was also calculated for reference. 

The cr(u)) spectra were derived from a function as 
follow in which direct interband transitions were as- 
sumed; 18 ) 
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Here, the \n'k) and \nk) states denote the unoccupied 
and occupied states, respectively, e and p are the po- 
larization of light and the momentum of the electron, 
respectively, /(e n s) is the Fermi-Dirac distribution func- 
tion, huj nn i = e n £ — e ,jg is the energy difference of the 
unoccupied and occupied states and T is the lifetime pa- 
rameter. In the calculation, T = 1 meV was assumed. 

We measured the reflectivity [-R(w)] spectra of CeSn 3 
and Celn3. The near- normal incident optical R(u>) spec- 
tra of polycrystalline CeSn3 and Celn3 were acquired 
in the photon energy rage of 2 meV - 30 eV using 
synchrotron-based equipment at UVSOR-II. The u(ui) 
spectra were derived from Kramers-Kronig analysis of 
the R(ui) spectrum. For CePd3, we used the a(w) spec- 
trum reported elsewhere. 13 ' 

The calculated band structures of CePd 3 with and 
without SOI are plotted in Fig. 1. The highly disper- 
sive bands below -1 eV, highly dispersive bands above 
-1 eV, and flat bands at - 0.5 eV above E F mainly 
originate from the Pd 4d, nearly free electron, and Ce 4/ 
states, respectively. In Fig. 1(a), the Ce 4/ state splits 
into two bands with energy splitting of about 0.3 eV by 
SOI. The band structure is fundamentally identical to 
that previously reported. 20 ' 21 ' 

The <t(oj) spectra calculated from the two band struc- 
tures with and without SOI in Fig. 1 are plotted in Fig. 2. 
The experimental a(u>) spectrum of CePd3 at T = 8 K 
is also plotted in the figure. The experimental spectrum, 
has three large peaks at 0.25, 0.55, and 0.78 eV, and one 
small peak at 0.04 eV. The calculated ct(uj) spectrum 
with SOI also has three peaks at 0.3, 0.56, and 0.75 eV 
and one shoulder at 0.1 eV, despite the very weak a(u>) 
intensity of LaPd3. This structure can be assigned to 
the experimental structure as follows: The peaks mainly 
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Fig. 2. (Color online) Experimental and calculated optical con- 
ductivity [<x(o>)] spectra of CePd3. The calculated ct(lu) spectrum 
of LaPd3 is also plotted for reference. Two different calculations 
with and without SOI were performed. The experimental spec- 
trum, which was obtained at T = 8 K, is derived from that 
reported elsewhere. 13 ' 



appear in the T-X-R-M plane and the initial state of 
all peaks is commonly the highly dispersive conduction 
band (the line with solid circles in Fig. 1). The peaks at 
0.3 eV and 0.78 eV in the calculated spectrum originate 
from the optical transitions to the Ce 4/7 / 2 state near the 
T-point and near the bottom of the A-axis, respectively. 
The small peak at 0.55 eV originates from the transition 
near the A- and S-axes. The shoulder structure at 0.1 eV 
in the calculation originates from the transition from the 
occupied flat band to the unoccupied Ce 4/ 5 / 2 state near 
the T-point. Therefore, all of the mid-IR peaks of CePd3 
can be explained by the band structure calculation with 
SOI. 

In the calculated spectrum without SOI, one large 
peak at 0.7 eV appears, as shown in Fig. 2. The cr(u>) 
spectrum is fundamentally equal to the calculation re- 
ported previously. 17 ' The spectrum cannot reproduce the 
experimental triple-peak structure. The SOS of the Ce 4/ 
state therefore plays an important role in the mid-IR 
triple-peak structure of CePds; i.e., the mid-IR peaks 
originate from the band structure including the SOS of 
the Ce 4/ state. 

To investigate the origin of the mid-IR peak in other 
Ce compounds, the band structure of CeSn3, a mixed va- 
lent material, was calculated as shown in Fig. 3 and the 
calculated cr(u>) spectrum was compared with the experi- 
mental spectrum, as shown in Fig. 4. The band structure 
is consistent with that previously reported by Hasegawa 
et al. 22 } The experimental spectrum has two peaks at 
about 0.15 and 0.45 eV, and a large Drude structure 
due to the carriers. Since the calculated c(w) spectrum 
only has the interband component, the experimental in- 
terband spectrum after subtracting the Drude part from 
the experimental cr(u>) spectrum is compared with the 
calculation. The Drude part was evaluated by the slope 
of the experimental spectrum at around 0.04 eV (see the 
inset of Fig. 4). The calculated <j(lo) spectrum has a large 
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Fig. 3. Calculated electronic structure of CeSri3 with the spin- 
orbit interaction. The bands with solid circles, open triangles, 
solid squares and open circles are the origin of the mid-IR peaks 
in this material. 



Fig. 5. Calculated electronic structure of Celn3 with the spin- 
orbit interaction. The bands with bold lines relates to the mid-IR 
peak of Celn3 . 




Fig. 4. (Color online) Experimental and calculated optical con- 
ductivity [o"((^)] spectra of CeSn3. The calculation was performed 
with SOI. The experimental spectrum was recorded at T = 8 K. 
The interband spectrum was obtained by subtracting the Drude 
component due to the carriers from the experimental a-(ui) spec- 
trum. The inset shows the lower energy part of the experimental 
(t(oj) spectrum and the Drude curve fitted at the slope at around 
0.04 eV. 



peak at 0.15 eV, a dip at 0.4 eV, and a shoulder at 0.5 eV, 
and is in good agreement with the experimental spec- 
trum. According to the calculation for each band contri- 
bution (not shown), the peak at 0.15 eV can be assigned 
to originate from the transition between the similarly dis- 
persive bands of solid circle and open triangle in Fig. 3 in 
between the M (101) and M (HO)-points. The two bands 
expand with the parallel dispersion in a wide area near 
the M-point, so a large joint density of states or a a(oj) 
peak is produced. On the other hand, the shoulder struc- 
ture at 0.5 eV originates from the transition between the 
bands indicated by solid circles and open circles in Fig. 3. 
In addition, small peaks at around 0.2 eV, which origi- 
nates from the transition between the bands of solid cir- 
cles and solid squares, also correspond to the shoulder at 
0.2 eV in the experimental spectrum. The experimental 



mid-IR peak structure of CeSn3 can therefore also be ex- 
plained by the band structure calculation. However, it is 
noted that the intensity of the calculated ct(cj) spectrum 
is more than double the experimental value. The same 
tendency has been observed in other materials, although 
the reason for this is not clear at present. 23 ) One possibil- 
ity is the strong electron correlation in Ce 4/ electrons. 
Since the value of the on-site Coulomb repulsion energy 
(U) is about 6 eV in Ce 4/ electrons, if the cf hybridiza- 
tion intensity is not so large, some of the Ce 4/ states 
move away from E-p. 24 ^ Therefore, the transition inten- 
sity to the Ce 4/ state near Ep may become smaller that 
the calculated value. To confirm the electron correlation 
effect of the mid-IR peak, band calculation incorporating 
dynamical mean-field theory is needed. 

Finally, we investigate whether the mid-IR peak of 
Celn3 in the antiferromagnetic ground state can be ex- 
plained by the band calculation. Generally, the band 
structure calculation of a material with local 4/ elec- 
trons is performed with U, namely LDA+U method. 18 ) 
In the case of Celn3, however, the U is considered to 
be small because the material is located near QCP. 
Actually, the angle-resolved photoemission spectra of 
CeRu2(Sio.82Geo.i8)2 that is located in the slightly lo- 
cal regime near QCP can be explained by the itinerant 
band calculation. 25 ) Since Celri3 is considered to be in 
the same situation as CeRu2(Sin.82Geo.i8)2, the itinerant 
band calculation of Celn3 was performed and compared 
with the experimental cr(uj) spectrum. 

The calculated electronic structure shown in Fig. 5 
is consistent with the previous work. 26 ) Fi gure 6 shows 
the calculated and experimental <t(cj) spectra of Celn3. 
To make a comparison with the paramagnetic band cal- 
culation, the experimental spectrum was taken in the 
paramagnetic phase at T — 60 K. The interband part 
after subtracting the Drude curve, which was fitted at 
around 0.04 eV, from the experimental <j{uS) spectrum 
has a broad mid-IR peak at 0.2 eV. The peak does not 
have the SOS structure, because the peak is broader 
than the splitting energy. The calculated spectrum has 
some small peaks originating from the optical transition 
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Fig. 6. (Color online) Experimental and calculated optical con- 
ductivity [o-(cj)] spectra of Celn3 . The calculation was performed 
with SOI. The experimental spectrum was taken in the paramag- 
netic phase at T = 60 K. The interband spectrum was obtained 
by subtracting the Drude component due to the carriers from the 
experimental o(lo) spectrum. The Drude curve was evaluated at 
the slope at around 0.04 eV. The inset shows the calculated o(u) 
spectrum compared with the experimental interband spectrum 
in the wide energy range below 4 eV. 



in some bands (bold lines in Fig. 5) near Ep. Although 
the calculated spectrum has some small peaks, the over- 
all spectral shape corresponds to the experimental spec- 
trum below 0.8 eV, for instance the flat-top shape of 
the experimental spectrum from 0.1 to 0.4 eV. Mean- 
while, several peaks in the energy range of 0.9 - 1.9 eV 
in the calculation can be regarded to be equivalent to 
the broad peak at 1.4 eV in the experiment as shown 
in the inset of Fig. 6, because the experimental life time 
broadening seems to be larger than that in the calculated 
spectrum in this energy region. So the dip structure at 
around 0.6 eV in the calculated c(w) spectrum also cor- 
responds to the minimum at 0.9 eV in the experimental 
spectrum. The broad mid-IR peak of Celn3 can therefore 
also be concluded to originate from the band structure. 
As mentioned earlier, the experimental spectrum in the 
antifcrromagnetic phase at 8 K is identical to that in the 
paramagnetic phase at 60 K. This result is also evidence 
that the mid-IR peak originates from the band structure. 

So far, the experimental <j{tu) spectra of CePd3, CeSn 3 , 
and Celn3 have been demonstrated to be in good agree- 
ment with the spectra derived from the band structure 
calculation including SOI. In addition, the triple-peak 
structure in CePd3 can be explained with the SOS of the 
Ce 4/ state. These results indicate that the mid-IR peaks 
of these materials fundamentally originate from the band 
structure and the SOS of the Ce 4/ state plays an im- 
portant role in the shape of the mid-IR peaks. This is 
consistent with the fact that the mid-IR peaks of many 
other Ce compounds have a double-peak structure with 
energy splitting of 0.25 eV. The spectral shapes of CePd3, 
CeSn3, and Celn3 are significantly different from one an- 
other. This indicates that the spectral shape originates 
from the particular electronic structure near Ep in each 
case. Therefore, the mid-IR peaks can be explained by 



the band structure calculation with SOS and without ref- 
erence to the cf hybridization gap. This result is reason- 
able in the mixed-valent Ce compounds, because the 4/ 
states can be treated as itinerant bands, for instance, the 
band structure calculation of CeSn3 can reproduce the 
angle-dependent frequency of the de Haas-van Alphcn 
effect. 22 ) 

To summarize, we calculated the <t(uj) spectra of 
CePd3, CeSn3 and Celn3 in the infrared region from 
the first-principle band calculation with SOI and com- 
pared them with the experimental spectra. The calcu- 
lated <j(uj) spectra are consistent with the experimental 
spectra, even though the experimental u{lo) spectra of 
these materials are significantly different from one an- 
other. The different spectral shapes originate from the 
different electronic structures near E-p. These results in- 
dicate that the mid-IR peaks originate not from the con- 
ventional cf hybridization gap but from the characteris- 
tic band structure. 
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